CERN-TH/2003-222 1 



m 
o 
o 

(N 

CD 

(N 
(N 



> 

O 
m 
o 



X 



QCD at zero baryon density 
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While the grand canonical partition function ZGcitA with chemical potential ji explicitly breaks the sym- 
metry with the Dirac determinant, the canonical partition function at fixed baryon number Zc{B) is manifestly 
Za-symmetric. We compare ZacilJ, = 0) and Zc{B = 0) formally and by numerical simulations, in particular 
with respect to properties of the deconfinement transition. Differences between the two ensembles, for physical 
observables characterising the phase transition, vanish with increasing lattice size. We show numerically that the 
free energy density is the same for both ensembles in the thermodynamic limit. 
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1. INTRODUCTION 

The grand canonical partition function with re- 
spect to the quark number 



ZGc{^J)^ / [i?C/]e-^«[^ldetM(C/;/^), 



(1) 



which is commonly used, explicitly breaks the Z3- 
centre symmetry. We compare it with the canon- 
ical partition function at fixed baryon number B 
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Figure 1. Distribution of the complex Polyakov 
loop in the grand canonical (top) and canonical 



Zc{B) = — dflie''^^'"ZGc{tJ' = ifJ-iT), (2) (bottom) ensembles, left: 4^ x 4, right: 6 
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which preserves the Z^-center symmetry. Here, 
we consider the case n = B = only. While 
in one ensemble the expectation value of the 
Polyakov loop is non-zero at any temperature, in 
the other it is always zero, as visible in Fig. 
Therefore, an important question arises: Do both 
formulations agree in the thermodynamic limit 
for physically relevant observables? This is a sub- 
tle issue, leading to a controversy about contribu- 
tions of non-zero triality^ states. Refs.^ suggest 
that non-zero triality ( "unphysical" , fractional B) 
states give a non-vanishing contribution even in 
the thermodynamic limit, and must be taken into 



*Talk presented by S. Kratochvila 

^The difference of the quark and anti-quark numbers mod- 
ulo 3, t = (#Q — #Q) mod 3, is called the quark triality. 



4. In the thermodynamic limit, the distributions 
agree for both ensembles, up to two additional 
Zs-rotations in the canonical ensemble. 

account. In contrast, f21 shows, for QCD with 
staggered fermions, that the only relevant sector 
of the grand canonical /i = description in the 
thermodynamic limit is the B = sector (the free 
energy densities of both descriptions are identi- 
cal for V — > 00). In other words, the canonical 
ensemble Zc{B = 0), which contains only zero- 
triality states, describes the same physics in the 
thermodynamic limit as the grand canonical en- 
semble ZQcifJ- = 0)7 which has contributions from 
all triality-sectors. Therefore, non-zero triality 
states do not contribute in this limit. By an ex- 
plicit comparison of Zqc{^i = 0) with Zc{B = 0), 
we hope to settle this discussion. 
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Figure 2. Phase diagram of Zcciii^i) in the 
(/i/, T)-plane. The transitions (fuU hnes) are first 
order for four flavours of staggered quarks with 
mass ma — 0.05. 

2. THE CANONICAL ENSEMBLE 

We briefly recall the construction of the canon- 
ical partition function Zc{B). First, one fixes the 
number of quarks N = J (Px 4'{x) 70 i^i^) to Q 
by inserting a i5-function in the grand canonical 
partition function: 

Zc[Q) = I [DU][d^][d^]e-^^^^^-^^^^''^'^^5[N ~Q) 



The (5-function admits a Fourier representation: 



Zc{Q) = ^ / 



-iQfj.1 



with the new variable flj. One recognizes /i/ — 
jliT as an imaginary chemical potential, so that 



ZciQ) = 



1 

2^ 



■TZaciif^i)- 



The grand canonical partition function Zcci^fJ-i) 
as a function of an imaginary chemical potential 
has interesting properties [H], see Fig. [3 It is even 
and 27rT/3-periodic in /j,/. For our purposes, the 
most important consequence is, that Zc{Q ^ 
mod 3) = 0. Therefore, we define the baryon 
number B = 3Q and end up with Eq.Q. 

3. RESULTS 

We simulate QCD with Nf — 4 flavours of stag- 
gered quarks of mass ma ~ 0.05, using Hybrid 
Monte Carlo. We measure the Polyakov loop, the 



plaquette and the chiral condensate on lattices 
with extents 4^, G'^ and 8"^ x 4 at seven different 
/3's, from f3i — 4.9 to — 5.1. We analyse the 
results using multihistogram reweighting W^. 

In order to sample the B = canon- 
ical ensemble, we supplement the ordinary 
Hybrid Monte Carlo at fixed /i/ with a 
Metropolis update fij —>■ fx'j, with acceptance 

min ( 1, ■ 



dot f { ip{fj.j)+m) ^ rp-j^g determinant ratio 



dot" f{I?{lij)+m), 

is evaluated with a stochastic estimator, namely 
(e-K mp'i)+^r"f^^( i/>M+^)"f^^v\^+\n\^)^ where 
?7 is a gaussian complex vector. In addi- 
tion, one can perform a "Z3 move" at any 



time (with 
Ui{x,t = 0) 



acceptance 1): fJ^i — 
U4{x,t = 0)e=F^if , Vx. 
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Figure 3. Susceptibility of ';/'■!/' on all analysed 
lattices and both ensembles, left: Zccif^ = 0), 
right: Zc{B — 0). Even for the smallest, 4**, 
lattice, differences are barely visible. 



Fig. 3 shows the susceptibility of the chiral con- 
densate ipip in both ensembles. The peak shifts 
slightly in f3 for the 4^^ x 4 lattices. For larger 
lattice sizes, the difference disappears and both 
ensembles indicate the same critical (3c already 
for the 8"^ x 4-lattices. We observe the same be- 
haviour for the average plaquette. This is our first 
evidence to answer the question about the infiu- 
ence of non-zero triality states: In the thermo- 
dynamic limit, non-zero triality sectors are sup- 
pressed and do not affect the results. The effect 
of non-zero triality sectors is smaller than the sta- 
tistical errors. 

In Fig. 4 we show a finite size analysis of observ- 
ables in the two ensembles. We plot the minimum 
of the Binder cumulant Cb{0) = 1 — as 
a function of the inverse volume \/V. For both 
the average plaquette and the chiral condensate. 
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Figure 4. Binder cumulant minimum versus in- 
verse volume for both ensembles (slightly shifted), 
left: average plaquette, right: chiral condensate. 
The thermodynamic extrapolation does not tend 
to |, indicating a first order transition. 



the thermodynamic extrapolation does not tend 
to I (top-limit of the plot) - indicative of a first 
order phase transition (as known from the litera- 
ture [5]). Furthermore, the slopes of the linear fit 
are identical in both ensembles within statistical 
errors. Therefore, we cannot even claim smaller 
finite size effects in one of the two ensembles. 
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Figure 5. Free energy density as a function of /i/. 
left side: (3 = 4.95 < /3c, right: /3 = 5.10 > /3c. 
In the latter case, a parabola nicely describes the 
data for the larger lattices. 

In Fig. 5 we show the free energy density versus 
/i/ for both /3 < /3c and /3 > /3c. In both cases, 
we clearly observe a minimum at /i/ = 0. There- 
fore, in the thermodynamic limit, only — 
mod will survive. This establishes the equiv- 
alence of Zc{B = 0) with ZGci^J■ = 0). For 
/3 = 4.95 < /3c, we observe that the free energy 
density at /u/ = zLttT/ 3 seems to have zero deriva- 
tive, indicating a crossover (as expected from the 
phase diagram Fig. [SJ. Instead, for /3 = 5.10 > /3c 
we expect a cusp to develop, due to the first order 
phase transition. Indeed, it appears likely as the 
volume increases. A simple parabola gives a good 
description of the data, with a curvature in line 
with expectations from high temperature [Hj. 



4. CONCLUSION 

The thermodynamic limits of Zqc{pl = 0) and 
Zc{B — 0) seem to agree with respect to the de- 
confinement phase transition. The effect of non- 
zero triality sectors is smaller than the statisti- 
cal errors, already for S'^ x 4 lattices. Therefore, 
these states can be included or excluded without 
affecting the results in the thermodynamic limit. 
Moreover, finite size effects are equivalent even on 
the smallest, 4"^ x 4, lattice. 

The canonical formulation requires a centre- 
symmetric formulation of QCD, which can be 
achieved very simply with negligible computer 
overhead, by adding a single d.o.f. /u/ up- 
dated by Metropolis. Hence, the presence of a 
fermionic determinant does not necessarily break 
the centre-symmetry. 

Note, that another grand canonical partition 
function Zqc{^) can be built from the Zc(_B)'s, 
which is completely centre-symmetric. This par- 
tition function will give identical expectation 
values to the usual ZacifJ-), except for centre- 
sensitive observables (Polyakov loop), whose dif- 
ferent behaviour can be understood by invoking 
spontaneous center-symmetry breaking \7\. 
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